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The central engines of two unusual radio- intermediate/quiet 
active galactic nuclei: III Zw 2 and PG 1407+265 

Liang Chen 1 , Xinwu Cao 1 and J. M. Bai 2,3 
ABSTRACT 

We use the accretion disk/corona+jet model to fit the multi-band spectral 
energy distributions (SEDs) of two unusual radio- intermediate/quiet quasars. It 
is found that the optical/UV emission of III Zw 2 is probably dominated by the 
emission from the accretion disk. The X-ray emission should be dominated by 
the radiation from the jet, while the contribution of the disk corona is negligible. 
The optical/UV component in the SED of PG 1407+265 can be well modeled 
as the emission from the accretion disk, while the IR component is attributed to 
the thermal radiation from the dust torus with an opening angle ~ 50°. If the 
X-ray continuum emission is dominated by the synchrotron emission of the jet, 
the source should be a "high peak frequency blazar", which obviously deviates 
the normal blazar sequence. The observed SED can also be fitted quite well 
by the accretion disk/corona model with the viscosity parameter a = 0.5. The 
spectrum of the accretion disk/corona in PG 1407+265 satisfies the weak line 
quasar criterion suggested in Laor & Davis. 
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Introduction 



In optically selected quasar samples, quasars with similar optical properties exhibit very 
different properties in radio bands. The radio loudness parameter R, the ratio of the radio 
flux at 5 GHz to the optical fl ux at B-band (R = fr,c M r J fn), is used as an ind icator of 



the radio properties of quasars (iKellermann et al.lll989l ). iKellermann et al.l (119891 ) found a 



dichotomy in the radio loudness distribution for an optically selected quasar sample. The 
quasars with radio loudness R > 10 are defined as radio-loud (RL) quasars, while those with 
R < 10 are named as radio-quiet (RQ) quasars. It was found that RL active galactic nuclei 
(AGNs) have powerful relativistic jets, while most of RQ AGNs have no or very weak jets. 
The distribution of radio-loudness may provide useful clues on the jet formation mechanism 
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). The 



FIRST (Faint Images of the Radio Sky at Twenty Centimeters) detected quasars showed tha t 
the radio loudness distribution is not bimodal, but rather continuous fjWhite et al.l 120001 ) . 
It is still debating on whether the radio -loudness distribution is bimodal or continuous (see, 
Ivezic et al.lbooi Icirasuolo et al.lbooih . 



The radio-loudness p arameter can be a good indicator of t he jet properties of AGNs 



in the sense of statistics. iMiller et al.l (119931 ) and iFalcke et al.l (119961 ) identified a number 



of quasars with "radio intermediate loudness". They suggested that these sources might 
be relativistically boosted radio- weak quasars ("radio-weak blazars"). In fact, the very long 



AGNs exhibit jet structure and high-brightness temperature radio cores 
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20(3). In this work, we choose two unusual RQ/radio intermediate sources with 



relativistic jets, III Zw 2 and PG 1407+265, to explore the physics of the disk-jet connection 
in this kind of the sources. As these two sources have both RL and RQ characteristics, 
the present investigation may reveal the relation of the disk-jet connection with the radio- 
loudness parameter, which may provide clues on the origin of the radio-loudness distribution. 



Arpl (119681 ) classified III Zw 2 (PG 0007+106, Mrk 1501, z = 0.089) as a Seyfert I galaxy. 



which is also i ncluded in the PG quasar sam ple (jSchmidt fc Greenlll983l ). It is hosted in a 
spiral galaxy (IHutchings &: Campbelll Il983l ). which is typical characteristics of RQ AGNs. 
The extend ed radio emission i s very weak compared w ith its core emission (~ 50 — 100 mJy 
at 1.4 GHz, lUnger et al.lll987l ; iBrunthaler et al.ll2005l ). and superluminal motions of the jet 
compon ents in this source hay e been detected in the VLBA (Very Long Baseline Array) 
images (IBrunthaler et al.ll2000l ). The apparent velocity of the moving jet components is 
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/3app = 1-25 ± 0.09 in units of light speed is derived from the VLBA images at 43 GHz. Ill 
Zw 2 exhibits violent variability throughout a ll wavebands (for e xample, more than an order 
of magnitude of variability in radio bands, see lSalvi et al.ll2002l ). IChen et al.l ( 120101 ) modeled 
the multi-waveband SED of this sources, and suggested that III Zw 2 is a possible 7-ray 
source and could be detected by the Fermi /LAT in the future. This source pos sesses typical 



blazar-like properties, however, it is included in the radio-intermediate sample (IFalcke et al. 



1996|). 



PG 1407+265 is a RQ AGN (z = 0.94, Kellermann et al.lll989l : McDowell et al.lll995l ) 
with some unusual properties. Its flux ratios in the r adio, X-ray, and optica l wavebands are 
typical of normal RQ quasars (IPlotkin et al.ll2010l ). iMcDowell et al.l ( 119951 ) found that the 
emission lines of PG 1407+265 have very small equivalent widths (except for Ha), but the 
full width at half- maximum (FWHM) reaches v ~ 7000 — 12000 km s _1 . This is a weak 
line quasar ( WLQ). The nature of WLQs is still unclear so far, though they are extensively 



studied (see. lLaor fc David 1201 ll . and the references therein). The optical-to-X-ray spectral 



index, a ox , of ty pical RQ AGNs rang es from 1.2 to 1.8, and the slope steepens with increasing 
UV luminosity ( jStrateva et al.ll2005l ). During the high state of PG 1407+265, the optical-to- 
X-ray spectrum is flat (a ox =1.09), and there are larger fluctuatio ns in the X-r ay bands than 
the UV bands (i.e., the slope flattens with increasing luminosity, lGalloll2006l ). which seems 
to be different from that in typical RQ AGNs. It is fou nd that the UV variability correlates 
with that of the X-ray emission in the high state, and iGallol (120061 ) suggested that the soft 
X-rays may be the jet's non-thermal emission in the high state. This implies that both the 
emission in the optical and X-ray bands may come from the jets. It is interesting to note 
that the source PG 1407+265 may probably have a relativistic jet moving toward us with a 
Doppler factor of S > 10 from its radio morphology and the variations in radio wavebands 
felundell et aDl2003[ ). To our knowle dge, PG 1407+265 i s the only RQ AGN with relativistic 
jets reported in the literature so far. I Chen fc Bail (120111 ) fitted the multi-waveband SED of 
this source with the homogeneous sphere jet model, and they found that the jet is required 
to be in very extreme conditions. Thus, they suggested that the radiation of the accretion 
disk may not be negligible in the optical/UV bands of the SED for this source. 

In this paper, we adopt an accretion disk/corona+jet model to explore the disk-jet 
connection in these two unusual AGNs. In Section 2, we give a brief summary of the jet 
and the disk-corona emission models used to fit the observed SEDs. Section 3 contains the 
results. The discussion and conclusions are given in Sections 4 and 5. Throughout this 
paper, a cosmology with H Q = 70 km s~ 1 Mpc~ 1 , Q m = 0.3 and f2 A = 0.7 is adopted. 
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The model 



The homogeneous sphere jet model was widely used to explain the observed SEDs 
of bl azars, with which the main f eatures of the blazars' sp ectra can be well reproduced 



(e.g., 



Tavecchio et al.l Il998l . l200ll ; ICelotti fc Ghisellinil 120081 ). For RQ AGNs, the accre- 
tion disk/corona model was d eveloped to explain the optical/UV and X-ray spectra (e.g., 



Haardt fc Maraschilll99ll . ll993l ). In this work, we adopt the homogeneous sphere jet+ accretion 
disk/corona model to fit the observed SEDs of these two AGNs, and their physical properties 
can be derived. We briefly describe the jet, and the accretion disk/corona models employed 
in this work as follows. 



2.1. Jet model 



In this paper, we use one zone synchrotron+inverse Comptonization (IC) models to 
calculate the jet emission of HI Zw 2 and PG 1 407+265. The model was widely used in 
blazar SED modeling (e.g.. iGhisellini et al.ll2010l . and the references therein). The emission 
region i s assumed to be a homogeneous sphere with radius R embedded in the magne tic field 



Tavecchio et al. 1998, 2001; Celotti fc Ghisellini 2008; Ghisellini et al. 2010, for the 



B (see, 

details). A broken power law electron energy distribution 



N(i) 



7min < 7 < 70 

7o < 7 < 7 max; 



1962; Sikora et al. 1994; Inoue & Takahara 1996; Kirk et al. 1998; Ghisellini et all 



the blob (see, 


Kardashev 


Ghisellini et al. 


19981. for 



the detailed discussion). The parameters of this model include, the radius R of the blob, 
the magnetic field strength B, electron break energy 70, the minimum and maximum energy 
7min, 7max, of the electrons, the normalization of the particle number density N , and the 
indexes pi^ of the broken power law particle distribution. The observed spectrum of the jet 
can be calculated when the Lorentz factor T — 1/ \/l — ft 2 , the viewing angle 6 of the jet with 
respect with the line of sight, and the spectrum of the external seed photons, are supplied. 
The frequency and luminosity can be transformed from the jet frame to observational frame 
as: v = bv' I (1 + z) and uL u = 5 A i , 'L' l/ , 1 where the Doppler factor 5 = 1/ [r (1 — ft cos 9)], and 
the prime represents the value measured in the jet frame. The synchrotron self-absorption 
and the Klein-Nishin a effect in the inverse Compton scattering are proper ly considered in 
our calculations (see, iRybicki fe Lightmanlll979l ; iBlumenthal fc Gouldlll970l ). Both the self- 
synchrotron Compton (SSC) scattering and external Compton (EC) scattering are included 
in the calculation of the Compton scattering in the blob. 
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In this homogeneous sphere model, the jet power can be calculated if all the physical 
quantities of the sphere are specified, 



ttR pT cU{- ot , 

where the total energy density measured in the rest frame of the blob, 



(2) 



UL 



The energy density for electrons U' e = m e c 2 J ^(7)7^7, while the proton en ergy density U' r = 
Uliny - n/rrip) / (7) if charge neutra lity for pure hydrogen plasma is assumed (jCelotti fc Fabian 



1993 



Celotti fc Ghisellinil 120081 ). The broad waveband SED from high freque ncy radio emis- 



sion to 7-ray bands can be m odeled with this homogeneous sphere model (e.g.. lCelotti fc Ghisellini 



20081 ; iGhisellini et al.ll2010l ). The emission from blazars in low frequen cy radio band may 



dominantly be rad i ated from the inner conical jet near the b lack hole (e.g.. lBlandford fc Konigl 
19791 ; lKonigllll98ll ; IGhisellini et al.lll985l ; iJiang et al.lll998l ). which is beyond the scope of this 
work. In this work, we use the conventional homogeneous sphere model to fit the observed 
SEDs from high fre quency radio emission to 7-ray bands, which is similar to most of the 



previous works (e.g., ICelotti fc Ghisellinil l2008t IGhisellini et al.ll2010l ). 



2.2. Accretion disk/corona model 



The observed UV/optical continuum emission of AGNs is thought to be the ther- 
mal emission from the standard geometrically thin, optic ally thick accretion disks (e.g., 
Shields! 1 19781 ; iMalkan &; Sargentlll982l ; ISun &; Malkanlll989l ). while the observed power-law 
hard X-ray spectra of RQ AGNs are most likely due to the inverse-Compton scattering of 



soft photons on a population of hot electrons in the corona above the disk (IGaleev et al. 



19791 ; lHaardt fc Maraschilll99ll . |l993j). In the accretion disk-corona model, such soft pho- 
tons are from the cold disk, a fraction of which are Compton scattered by the hot elec- 
trons in the corona above the cold disk to the hard X-ray energy band. The disk-corona 



model was extensively explored in many previous works fe.g..|Haardt fc Maraschilll99ll . 11993 
Svensson fc Zdziarskilfl9~9ll ; IrCawaguchi et al.lboOll : iLiu et al.lbooi ICaoll2009h . In this disk 



corona scenario, most of the gravitational energy of the accreting matter is released in the 
cold disk, and a fraction of which is transported into the corona probably by magnetic fields. 
The magnetic fields generated in the cold disk are strongly buoyant, and a substantial frac- 
tion of magnetic energy is trans ported vertically t o heat the corona above the disk with the 



reconnection of the fi elds (e.g., |Pi Matted Il998l ; |Pi Matteo et al.l Il999t Merloni fc Fabian 



200lL 120021 : ICaoll2009h . In this work, we adopt the model given in[Cao| fioS) to calculate 
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the spectrum of the accretion disk/corona system. We summarize the accretion disk-corona 
model in this sub-section. 

The gravitational power dissipated in unit surface area of the accretion disk is given by 



^dissi 



8^ 



MQ k (Rf 



1 - 



Ryn 



V2' 



(3) 



where M is the mass accretion rate of the disk, £\(-R) is the Keplerian velocity at radius R, 
and i?in = 3i?s (IShakura fc Sunyaevlll973l ). The Schwarzschild radius Rs = 2GM./c 2 , where 
M. is the black hole mass. The accretion disk luminosity is 



£disk = 4vr / Q+^RdR 



GM.M 

2R-, n 



(4) 



In the absence of the corona, the surface temperature T s of the disk as a function of radius R 
can be calculated with <rT s 4 = Qd issi , and the spectrum of the disk is ava ilable by integrating 
blackbody emissivity B V (R) over radius R (IShakura fc Sunyaevl Il973l ). In this case, the 
spectrum of the accretion disk can be derived when the black hole mass and the accretion 
rate are specified. 

In the accretion disk/corona system, the corona is assumed to be heated by the recon- 
nection of the magnetic fields generated by the buoyancy instability in the disk. The power 
dissipated in the corona is estimated with, 

B 2 



Q2 



8tt 



pi 



(5) 



where p m is the magnetic pressure in the disk, and v p is the velocity of the magnetic flux 
transported vertically in the disk ( iDi Matted Il998l ). The rising speed v p is assumed to be 
proportional to the internal Alfven velocity, i.e., v p = bv&, in which b is of the order of unity 
for extremely evacuated magnetic tubes. We adopt b = 1 in all our calculations of this work. 

The soft photons from the disk are Compton scattered by the hot electrons in the corona 
to X-ray bands, and about half of the scattered photons are intercepted by the disk. The 
reflection albedo a is relatively low, a ~ 0.1 — 0.2, a nd most of the inciden t photons from the 
corona are re-radiated as blackbody radiation (e.g.. IZdziarski et al.lll999l ). Thus, the energy 
equation for the cold disk is 



^dissi Q 



4^T d 4 isk 
3r ' 



(6) 



where Tdi S k is the effective temperature in the mid-plane of the disk, and r = r cs + Tg is the 
optical depth in vertical direction of the disk. In this work, we adopt a = 0.15 in all our 
calculations. 
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The detailed physics for generating magnetic fields in the accretion disk is still quite 
unclear, and there are three different magnetic stress tensors were usually adopted, 



Pn 



aptot 



(7) 



Shakura & Sunvaev 


1973; 


Sakimoto & Coroniti 


1981 



Taam fc Liru I1984J ) . The parameters of the disk-corona model include the black hole mass 
M m , the dimensionless mass accretion rate m (m = M/M-^dd, M^dd = -^Edd/?7radC 2 , and the 
convention al radiative efficiency r] Tau d = 0.1 is adopted), and the viscosity parameter a (see 
Caoll2009[ for the details). 



3. Results 



We search the literature for the observational data of the multi- waveband SEDs of these 



two sources from radio to X-ray bands (ISalvi et al.l |2002| ; iKaastra fc de Kortd Il988l . and 
#i?i>Q). Obviously, simultaneous broad band SEDs are desired to be used to model the disk- 
jet systems in these sources. For III Zw 2, we search the literature and collect a broad band 
SED, which is quasi-simultaneous within half a year (see Table HJ). This source shows long- 
term variability on time scale of a few years, which is longer than the temporal span of the 
quasi-simultaneous data. Therefore, we mainly adopt this SED for our model fitting. For the 
source PG 1407+265, we collect and list all X-ray observations in Table HJ We note that the 
GINGA observation (Jun. 1987) is roughly simultaneous with IR and optical observations 
(see Table [3] and H]) . Unfortunately, only 2-10 keV flux is available without spectral index. 
It can be seen that the 2-10 keV flux of observation on Jan. 17, 1981 by EINSTEIN is 
similar to that of GINGA observation. We therefore use EINSTEIN observation data in the 
quasi-simultaneous spectrum of this source (see Table HJ). 



3.1. Ill Zw 2 



We search the literature and collect the quasi- simultaneous SED of III Zw 2 from radio 
to X-r ay bands, which is plotted in Figure [j] (the data are listed in Table [2]). ISalvi et al. 
( 120021 ) estimated the central black hole mass M. = 1O 9 M from the width of the broad- 
line H/3 of this source. Superluminal motion of the jet component with apparent velocity 



1 http://ned. ipac.caltech.edu/ 
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0api> = 1-25 ± 0.09 in HI Zw 2 has been detected with VLBA observations (IBrunthaler et al. 
20001 ) . iPopovic et al.l (120031 ) derived the inclination angle of the emission line disk 9 = 12°±5° 



for III Zw 2 from its broad emission line profile. Assuming the jet to be perpendicular to the 
emission line disk, the Lorentz factor of the jet, T ~ 2.06, and the Doppler factor, 5 ~ 3.35, 
are derived. 

In the calculations of the spectra from the jets, we have considered both the SSC and EC 
mechanisms. The external soft seed photons from the broad line region (BLR) are considered 
in our calculations. The VLBA observatio ns and the variability in radio bands imply that the 
location of the blob is within the BLR (see Brunthaler et al. I hoooh . The radius of the emission 
regi on is estimated with t he minimum variabili t y tim e scale R ~ ct var 5/(l + z) ~ 1.5 x 10 16 
cm J Jang fc Miller! Il997h . iKaastra & de Kortel Jl988l ) estimated the size of the BLR « 10 18 
cm, and the photon energy density within the BLR U ext ~ 3.8 x 10 -4 erg cm -3 . Thus, the 
multi- waveband spectrum of the jet can be calculated if the values of the parameters, Nq, B, 
7o, 7min, Pi and P2, are provided. Besides the spectrum of the jet, we need to calculate the 
spectra from the accretion disk/corona system in this source. First, we consider the simplest 
case, i.e., a bare accretion disk without a corona. In this case, the spectrum of the accretion 
disk is only dependent on the mass accretion rate (see Section 2). 

From Figure [U it can be seen that the quasi-simultaneous SED seems to show two 
bumps at IR and UV bands respectively. The UV bump is a typical character of RQ AGN, 
which can be explained as the thermal emission from an accretion disk emission. We model 
the UV bump as disk emission, and jet emission accounts for the IR bump. In the model 
fitting the quasi-simultaneous SED, the least squares method is adopted to derive the best 
fitted model parameters. The model parameters, the accretion rate m, the minimum electron 
energy 7 m i n , and the magnetic field strength B, are carefully tuned in order to obtain the 
best model fitting to the quasi-simultaneous SED. The best fitting parameters are listed in 
Table [T] (see Figure [TJ. With the values of the parameters determined from the fitting of 
the SEDs, the jet power can be calculated with Equation ([2]). In this case, jet power is 



^sim ^ x ^q45 er g g i an( j j-jjg accre tion disk luminosity L s ^ k « 1.1 x 10 45 erg s 

We also repeat our calculations by adopting the accretion disk/corona model instead of 
the bare accretion disk model as described above. There are different choices of magnetic 
stresses (see Equation [7]) . Our results show that the contribution to the X-ray continuum 
spectra in this source is always negligible (see Tabled]). 

We plot the sen sitivity of the detector of Fermi /LAT in the Figure [TJ (purple lines, 
McEnery et al.ll2004l ). It can be seen that if the optical/UV emission is dominantly from the 
accretion disk, the predicted 7-ray flux is below the sensitivity of Fermi /LAT. 
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3.2. PG 1407+265 

The quasi- simult aneous SEP of PG 1407 +265 shows clear evidence of an optical-UV 



bump (see Figure [2]) . iHryniewicz et al.l ( 120101 ) estimated the black hole mass of this source 
to be 6 x 10 9 M Q . There are no (quasi-) simultaneous optical/IR spectra corresponding to 
the high state in the X-ray band. Thus, we have to lim it our model fittings to the quasi- 



simultaneous SED in the low state. iBlundell et al.l (120031 ) suggested that there is a relativistic 



jet in PG 1407+265, and the Doppler factor of the jet is estimated as 5 > 10. We adopt 
5 = 10 in our calculations for this source. The procedure of fitting its SED is the same as 
that for III Zw 2. We find that the optical/UV component in the SED may probably be 
from the accretion disk, while the origin of the IR component is still uncertain. There are 
two possibilities: the synchrotron emission from the jet, or the thermal radiation from the 
dust torus irradiated by the radiation of the central engine. In the case of the synchrotron 
emission, the minimum electron energy can be constrained, 7 m i n ~ 13.5, if equipartition 
between magnetic field and electron+proton energy is assumed. The jet power is Lj et w 
1.1 x 10 46 erg s _1 , while the disk luminosity is L^k ~ 6.8 x 10 46 erg s _1 . If the IR component 
is from the dust torus, our calculations require the minimum electron energy 7 m ; n 4.6, and 
the jet power is L- set rj 3.3 x 10 44 erg s" 1 , which is significantly lower than the disk luminosity. 
All the model parameters adopted for different cases are listed in Table |5j 

In order to explore the origin of the X-ray emission in PG 1407+265, we also model 
the SED of this source by including the contribution of the corona above the disk. In our 
calculations, we adopt different magnetic stresses respectively, and find that the contribution 
of the corona to the X-ray spectrum is always negligible, unless the stress r T(p = p m = ap tot 
is adopted with an relatively high viscosity parameter: a ~ 0.5 (see Figure [3]). 



3.3. Blazar sequence 



We compare these two sources with the well studied blazar sample ( ICelotti k, Ghisellini 



20081 ). In Figure HI the blazar sequence is plotted, and we find that these two sources roughly 



follow the correlation defined by the blazar sample, except the source PG 1407+265 if the 
IR component is from the dust torus and the X-ray emission is the synchrotron emission 
from the jet. We also compare the relation of jet power and disk lu minosity with tha t of th e 
blazar sample in Figure [5j in which the blazar sample is taken from iGhisellini et al.l ( 120101 ). 
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Discussion 



The simultaneous SEDs either in the high or low states are not available for III Zw 2, 
and only a quasi-simultaneous SED is obtained by searching the data in literature. This 
quasi-simultaneous SED shows a big blue bump, which is a typical character of RQ quasar 
and can be modeled quite well by the accretion disk model (see Figured]). Compared with 
RQ AGNs, it is known that relatively high polarization in the optical waveband is observed 
in blazars, because the emission in this waveband is dominated by the synchrotron emission 
from the relativistic jets. We note that almos t no polarization has been observed in the 
optical waveband for III Zw 2 (~ 0.28 ± 0.19%, iBerriman et al.lll990l). and the ampl itude of 
its optical variability is smaller than the radio and X-ray emission (ISalvi et al.ll2002l ). which 
imply that the optical/UV emission in this source may not be dominantly from the jet. We 
estimate the thermal timescale of the accretion d isk (assuming black hole spin a = 0, and 
using the same method as that in iLiu et al.l 120081 ). Thermal ~ 2 — 5 yr for this source in the 
low or the high states. This i s roughly co n sisten t with the observed optical/UV variability 
timescales (ISalvi et al.l 120021 ). IChen et al.l (120101 ) showed that the predicted 7-ray emission 
can be be detected by Fermi /LAT, if the jet emission is responsible for the optical/UV 
continuum emission. However, the model fittings to the simultaneous SED of this source in 
this work show that the UV/optical emission originates from the accretion disk (see Figure 
[TJ), and the 7-ray emission from the jet is below the sensitivity of Fermi / LAT. We note that 
III Zw 2 is not included in either 1LAC ( the first LAT AGN Ca talog, lAbdo et al.l boioh 
or 2LAC (the second LAT AGN Catalog, Ackermann et al.l boilh . We suggest that the 
optical/UV continuum spectra of this source may probably be dominantly emitted from the 
accretion disk. This issue can be sorted out if the simultaneous multi-waveband observations 
on this source both in low and high states are performed in the future. 

We also calculate the X-ray spectrum of the corona above the disk, and find that the 
contribution of the corona to the X-ray emission is always negligible, i. e., the X-ray c ontin- 
uum spectra are mainly emitted from the relativistic jet in this source. ISikora et al.l (120071 ) 
investigated the radio loudness of a sample with 199 sources consisting of broad-line radio 
galaxies (BLRGs), radio- loud quasars (RLQs), Seyfert galaxies (SGs), low-ionization nu- 
clear emission-line region galaxies (LINERs), and Fanaroff- Riley type I radio galaxies (FR I 
RGs). They found that there are two distinct, approximately parallel tracks in the plot of 
log R — log A (the Eddington ratio A = L bo i/L Edd , where L bo i is the bolometr ic luminosity 



and L khh is the Eddington luminosity). We compare III Zw 2 with the results in lSikora et al. 



( 120071 ). and find that this source locates between the RL and RQ trac ks (Eddington radio : 
A low « 0.019, A hig h ^ 0. 072, A sim « 0.039; radio-loudness R ~ 100 - 200. iSlcke et all f fl996h : 



Brunthaler et al.l fcOOOf l). 
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St rong Fe Ka line emission w ith equivalent width EW~ 220 A has been detected in III 
Zw 2 (jJimenez-Bailon et al.ll2005l ). It roughly follows the relation between the equivalent 
width and luminosity for RQ quasars, of whi ch the X-ray contin uum emission is supposed 
to be from the accretion disk/corona system (INandra et al.lll997l ). Our results suggest that 
only a small fraction of the observed X-ray continuum is from the corona of the accretion 
disk, which implies that this source in fact deviates significantly from the relation between 
the equivalent width and luminosity defined by the RQ quasar sample if the X-ray continuum 
emission from the jet is properly subtracted in this source. The X-ray emission from the 
corona is obviously not enough to power such a strong Fe Ka in this source, and therefore 
the detailed studies on X-ray reflection geometry in III Zw 2 may be necessary for resolving 
this issue, which is beyond the scope of this work. 

The SED of PG 1407+265 exhibits a clear component in optical/UV bands , which can 
be well fitted by the accretion disk model (see, Figure [2]). IChen fe Bail (120111 ) fitted the 
multi-band SED of this source using the one- zone jet model, and found that extreme values 
of some model parameters are required. Similar to III Zw 2, almost no polarizati on in optical 
emission has been detected in this source (~ 0.24 + 0. 16%. iBerriman et al.lll990l ). which may 
also suggest that most of the optical-UV emission may not come from the jet. 

The origin of the X-ray emission in this source is somewhat uncertain, which could 
be dominantly either from the jet or the corona above the disk. We firstly adopt the 
jet+accretion disk/corona model to fit its SEDs. In this case, our results show that the 
X-ray emission is dominantly from the jet. Unlike III Zw 2, a component in IR wave bands 
has been observed in PG 1407+265, which may originate from the synchrotron emission from 
the jet, or the thermal radiation from the dust torus irradiated by the radiation of the central 
engine. In the case of the synchrotron emission, the derived jet power is L )et w 1.1 x 10 46 
erg s -1 , which is lower than the disk luminosity (Ldisk ~ 6.8 x 10 46 erg s _1 ). If the IR com- 
ponent is alternatively assumed to be emitted from the dust torus, our calculations indicate 
that the jet power L- ]e ,. t ^ 3.3 x 1 44 erg s -1 , which is also significantly lower than the disk 
luminosity. iBlundell et al.l (120031 ) estimated the lower limit on the jet power in this source, 
~ 7 x 10 43 erg s _1 , by assuming minimum of the total energy density, in which the minimum 
magnetic field strength B min is adopted. Considering that the realistic field strength could 
be significantly higher than B min , their estimate of the lower limit on jet power is consistent 
with our results in the case of the IR emission being from the torus. The opening angle 
of the dust torus can be estimated from the ratio of IR to optical/UV fluxes by assuming 
that the radiat ion of the disk is absorbed by the dust torus and re-radiates in IR wavebands 



sec 



Cad 120051 . for the detailed discussion). The ratio of IR to optical/UV emission ~ 0.6 
for PG 1407+265, and the torus opening angle ~ 50 ° is inferred, which is typical for RQ 



AGNs (Elvis et al.lll994J : ICaol 120051 : IShang et al.l 1201 lh . We find that the X-ray emission is 
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dominated by the synchrotron emissions from the jet (see the lower panel of Figure [2]), and 
this source should be a "high peak frequency blazar". This source is luminous, and it is 
quite different compared with normal blazars, be cause the blazar sequence expects a high 



peak frequency bl azar when its luminosity is low ( IFossati et al.lll998 
Chen fc Bail 120 llh . 



Ghisellini et al. 1998 



In order to investigate the origin of the X-ray emission in this source, we alternatively 
adopt the accretion disk/corona model to fit the SED for this source. We find that the 
spectra in the optical/UV and X-ray bands can be fairly well fitted with the magnetic stress 
r Tip = aptot, provided m = 0.3 and a = 0.5 are adopted (see Figure [3]). However, we note 
that the variability of the X-ray emission in this source can be around an order of magnitude 
between its low and high states (see Figure[3]), which is more violent than normal RQ quasars. 
Future simultaneous multi-band (UV/optical+X-ray) observations on this source will help 
resolve this issue. In the accretion disk/corona model, the variability in the X-ray band may 
lead to variable optical-UV emission, and vice versa. The observational data show strong 
variability in the X-ray band of this source, while no evidence of similar variability is found 
in the optical/UV band. Further simultaneous multi-waveband observations is expected to 
attack th e nature of the X-ra y emission in this source. We compare PG 1407+265 with the 
results in ISikora et al.l (120071 ). and find that this source locates in the R Q track (Eddington 
radio A « 0.33; radio-loudness R « 0.44,3.43, iWilkes & Elvisl (119871 ) : Kellermann et al. 



( 119891 )). Therefore, PG 1407+265 could be a typical RQ quasar, but contains relativistic 
weak jets. 



The nature of WLQs is still unclear. lLaor fc Davisl ( 1201 ll ) proposed that the accretion 
disk temperature decreases with increasing black hole mass, and the fraction of the photons 
with the energy that can ionize the BLR decreases. This may account for the weak line 
emission in some AGNs. They suggested that, the quasar will be lineless if the fraction 
of the accretion disk emission above the frequency v = 3.29 x 10 15 Hz is less than 0.01 
(^>3.29xio 16 Hz/^ ^ 0.01), while a weak lin e quasar appears if £ „>3.29xio 15 Hz/£ < 0.1. PG 
1407+265 is the first WLQ studied in detail rtMcDowell et al.lll995h . We calculate the fraction 
of disk-corona emission above t he frequency, A/^2fl xio 15 Hz/-£ ~ 0.08 — 0.1, which satisfies 
the WLQ criterion suggested by lLaor &: Davisl (120111 ). 



We compare these two sources with the well studied blazar sample ICelotti fc Ghisellini 



( 120081 ) . We find that these two sources roughly follow the blazar sequence defined by the 
blazar sample, and also compare the relation of jet power and disk luminosity with that of 
blazars. We find that these two sources do not deviate much from the blazars, except the 
source PG 1407+265, if its IR component is assumed to be from the dust torus and the 
X-ray emission is the synchrotron emission from the jet. 
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5. Summary 

We summarize the main conclusions of this paper as follows. 

1. The optical/UV emission of III Zw 2 is dominantly from the accretion disk. The 
X-ray emission should be dominated by the radiation from the jet, while the contribution of 
the disk corona is negligible. 

2. The predicted 7-ray emission is below the sensitivity of Fermi /LAT, provided the 
UV is emission originated from the accretion disk. 

3. The optical/UV component in the SED of PG 1407+265 can be well modeled as 
the emission from the accretion disk, while the IR component is attributed to the thermal 
radiation from the dust torus with an opening angle ~ 50°. If the X-ray continuum emission 
is dominantly from the jet, the source should be a "high peak frequency blazar", which 
obviously deviates the normal blazar sequence. 

4. The observed SED can also be fitted quite well by the accretion disk/corona model 
with the viscosity parameter a = 0.5. The spectrum of the accretion disk /coron a in the 



WLQ PG 1407+265 satisfies the WLQ criterion suggested bv lLaor & Davisl (120111 ). 
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Fig. 1. — The quasi- simultaneous SED of III Zw 2 (see Table |2] for references of the data). 
The method of least squares is used to fit SED. The green solid line represents the emission 
from the jet, while the red solid line is for the disk/corona emission. The red dashed line 
represents the spectra from the bare accretion disk, while the green dashed line is for the 
emission from the jet in this case. The black line is the sum of the emission from the jet 
and the accr etion disk/corona. Th e purple line represents the sensitivity of 7-ray detector 
Fermi /LAT (IMcEnery et all 120041 ) . The parameters model calculations see Table [TJ 
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Fig. 2. — The same as Figured! but for the source of PG 1407+265. The method of least 
squares is used in fitting. 
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Fig. 3. — The same as Figure [2} but the SED is fitted with the accretion disk/corona model. 
The magnetic stress r Ttp = ap to t is used, and the best-fitted model parameters are: m = 0.3 
and a = 0.5. The solid line represents the spectra from the accretion disk/corona system. 
The dashed line represents the disk spectra, while the dotted line is for the spectra from the 
corona. 
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Fig. 4. — The blazar sequence : 7n v s. U( ot . The opened squares represent the blazars 
taken from ICelotti fc Ghisellinil (120081 ). The red stars represent III Zw 2, while the blue 
dots indicate PG 1407+265. The dot with a higher 7 corresponds to the case with the 
X-ray emission dominated by the synchrotron emission from the jet, while the lower one 
corresponds to the case of the X-ray emission emission dominated by the inverse Compton 
scattering in the jet. 
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Fig. 5.— Jet power Zvj c t vs. disk luminosity i^disk- The black triangles/squares present the 
blaz ars, while the triangles are for the sources only with upper limits on the disk luminosity 



see, 



Ghisellini et al.ll2010l ). The red stars represent III Zw 2, while the blue dots are for 



PG 1407+265. The dot with lower jet power corresponds to the case of the X-ray emission 
dominated by the synchrotron emission from the jet, while the higher one is for case of the 
X-ray emission dominated by the inverse Compton scattering in the jet. 
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Table 1. Parameters of III Zw 2 



Model 


B(Gs) 


N 


Pi 


P2 


7o 


Turin 


rh a 


T rV ^jet 


rt 

-^disk 


1 
2 


2.08 
2.08 


7.02x10 s 
7.02 x10 s 


2.3 

2.3 


7.0 
7.0 


1963 
1963 


20.0 
58.1 


5.1 x 10~ 3 ... 
1.0 x 10~ 2 0.2 


1.51 

a^/PgasPtot 0.40 


1.12 
1.25 



Note. — t £,j et anc j L disk are in unit of 10 45 erg/s. i? is the magnetic field, A*o the normalized 
number density, p\ 2 the indexes of the broken power law electron energy distribution, 70 the peak 
electron energy, 7 m ; n the minimum electron energies, m the dimensionless accretion rate, a the 
viscosity parameter. Model 1 indicates the dashed line in Figured] Model 2 is for the solid line in 
Figure [T] 



Table 2. quasi- simultaneous SED data of III Zw 



log^ log/„(Jy) Observation Date Ref. 



9.166 


-1.000±0.043 


Jan. 


1978 


S78 


9.689 


-0.569±0.016 


Jan. 


1978 


S78 


10.169 


0.114±0.033 


Jan. 


1978 


S78 


10.352 


0.415±0.033 


Jan. 


1978 


S78 


10.497 


0.663±0.038 


Jan. 


1978 


S78 


10.954 


0.839±0.044 


Jan. 


1978 


S78 


13.480 


-1.150±0.070 


Jul. 


1977 


N79 


13.934 


-1.450±0.060 


Jul. 


1977 


N79 


14.134 


-1.770±0.020 


Jul. 


1977 


N79 


14.260 


-2.010±0.020 


Jul. 


1977 


N79 


14.380 


-2.070±0.050 


Jul. 


1977 


N79 


14.450 


-2.120 


Jul. 


1977 


N79 


14.500 


-2.170 


Jul. 


1977 


N79 


14.550 


-2.230 


Jul. 


1977 


N79 


14.600 


-2.300 


Jul. 


1977 


N79 


14.650 


-2.350 


Jul. 


1977 


N79 


14.700 


-2.390 


Jul. 


1977 


N79 


14.750 


-2.430 


Jul. 


1977 


N79 


14.800 


-2.420 


Jul. 


1977 


N79 


14.850 


-2.430 


Jul. 


1977 


N79 


14.900 


-2.440 


Jul. 


1977 


N79 


14.950 


-2.450 


Jul. 


1977 


N79 


17.680 


-5.523 


Aug. 


1977 


K88 


18.379 


-5.732 


Aug. 


1977 


K88 



Note. — The quasi-simultaneous spectra of III 
Z w 2 as plotted in F i gureED The r eference S78 refers 



tolSchnopper et al.l (I1978T). N79: iNeueebauer et al 
(|l979h . K88: iKaastra fc de Kortel |l988h . 
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Table 3. SED data of PG 1407+265 



logi^ 


log/„(Jy) 


Observation Date 


Ref. 


9.17 


-2.06±0.03 


Sep. 


1993 


B96 


9.69 


-2.24±0.02 


Sep. 


1993 


B96 


9.93 


-2.33±0.02 


Sep. 


1993 


B96 


10.17 


-2.23±0.06 


Sep. 


1993 


B96 


12.69 


-0.77±0.13 


Jun. 


1996 


H03 


13.10 


-1.11±0.13 


Jun. 


1996 


H03 


13.61 


-1.77±0.13 


Jun. 


1996 


H03 


13.90 


-2.29±0.02 


Apr. 


1988 


E94,NED 


13.95 


-2.37±0.06 


Apr. 


1988 


E94,NED 


13.95 


-2.37±0.07 


Apr. 


1988 


E94,NED 


14.13 


-2.70±0.03 


Apr. 


1988 


E94,NED 


14.13 


-2.72±0.02 


Apr. 


1988 


E94,NED 


14.13 


-2.70±0.03 


Apr. 


1988 


E94,NED 


14.83 


-2.80±0.02 


May 


1986 


E94,NED 


14.73 


-2.73±0.01 


May 


1986 


E94,NED 


14.64 


-2.76±0.02 


May 


1986 


E94,NED 


14.53 


-2.78±0.02 


May 


1986 


E94,NED 


14.38 


-2.75±0.02 


May 


1986 


E94,NED 






May 


1986 


vqa wpn 

IL/c/^i, I V LjU 


14.26 


-2.79±0.02 


Apr. 


1988 


E94,NED 


14.26 


-2.75±0.02 


Apr. 


1988 


E94,NED 


17.38 


-6.36±0.17 






NED 


17.48 


-6.51 






NED 


17.50 


-5.47±0.01 






NED 


17.66 


-6.25±0.05 






NED 


17.99 


-7.00 






NED 


18.08 


-6.67 






NED 


18.10 


-6.80 






NED 


18.10 


-6.40 






NED 


18.16 


-7.26 






NED 



Note. 



SED 



data of 
and 



PG 



ence B96 refers to 


Barvainis et al. 


Haas et al. ( 


2003h. E94: Elvis et al. 


http://ned.ipac.caltech.edu/ 



1407+265 as 
The refer- 
dl996h . H03: 
dl994h. NED: 
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Table 4. X-ray spectra of PG 1407+265 



Observation Date fluxt Photon index /1-iokcV Telescope Ref. 



Jan. 17, 1981 


/ lkeV - O.Ut^l ^Jy 


r = 


2 2+ 1 ' 7 


1.27 


EINSTEIN 


W87 


Jun. 19-20, 1987 


/2-iokeV = 0.12 x 10- 11 erg/cm 2 /s 






1.2 


GINGA 


L97 


Jan. 19, 1992 


/ lkoV = l.oi&gj MJy 


r = 


2 61 +0 - 05 


1.62 


ROSAT 


M95 


Jul. 02, 1993 


/2-iokcV - l-38±g;g| x lO" 12 erg/cm 2 /s 


r = 


2 m+ 005 

^■ uo -0.05 


1.38 


ASCA 


GOO 


Jan. 23, 2001 


^2-iokcV = 8.85 x 10 45 erg/s 


r = 


qo+0.02 
^• oz -0.02 


1.33 


XMM-Newton 


G06 


Dec. 22, 2001 


/2-iokoV = 0.8 x 10~ 12 erg/cm 2 /s 


r = 


z - z ^-0.02 


0.8 


XMM-Newton 


F05 



Note. — X-ray spectra of PG 1407+265. t flux presented in references. * 2-10 ke V integral flux ( in uni t 
of 10 ~ 12 erg/cm 2 /s) ca l culate d from the flux* and the T. The re ference W87 r efers to I Wilkes fc Elvis (|1987l ) . 
L97: lLawson fc Turner! dl997t) . M95: iMcDowell et"ZI dl995h . GOO: iGeoree et all (l2000h . G06: iGallol (12001 . F05: 



Fang et al.l (|2005h . 



Table 5. Parameters of PG 1407+265 



Model 


B{Gs) 




w 


Pi 


P2 


7o 




m 


Q 






rt 

-^disk 


1 


29.4 


1.7xl0 5 


4.2 


1.8 


4.2 


95 


13.5 


0.08 


1.0 


a^PgasPtot 


11.3 


68.1 


2 


5.0 


1.8xl0 3 


4.2 


1.8 


4.2 


23148 


4.6 


0.08 


1.0 


ay/PgasPtot 


0.33 


68.1 


3 
















0.3 


0.5 


aptot 




144 



Note. — t R is the radius of emission region in unit of 10 15 cm. * Lj Ct and Ldisk are in unit 
of 10 45 erg/s. The model parameters have the same meanings as those in Table 1. Model 1 
corresponds to the upper panel of Figure O Model 2 is for the lower panel of Figure [21 Model 
corresponds to the result in Figure [3l 



